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Abstract 

We discuss the possible detection of a stochastic background of massive, non- 
relativistic scalar particles, through the cross correlation of the two LIGO 
interferometers in the initial, enhanced and advanced configuration. If the 
frequency corresponding to the mass of the scalar field lies in the detector 
sensitivity band, and the non-relativistic branch of the spectrum gives a sig- 
nificant contribution to energy density required to close the Universe, we find 
that the scalar background can induce a non-negligible signal, in competition 
with a possible signal produced by a stochastic background of gravitational 
radiation. 
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I. INTRODUCTION 



In the next few years many gravitational antennas will be collecting data, as the new 
interferometric detectors (GEO, LIGO, TAMA, VIRGO) |l|] will join the resonant detectors 
(ALLEGRO, AURIGA, EXPLORER, NAUTILUS, NIOBE) § already in operation, cov- 
ering a frequency band from ~ 10 Hz up to lKHz. Through the cross-correlation of these 
antennas it will be possible to search for a stochastic background of primordial gravitational 
radiation 0, [f|]. The detection of a relic background of gravitational waves would allow us 
to reconstruct the very early stages of cosmological evolution || and even an upper limit 
would be significant, since there are predictions of very high relic backgrounds || in principle 
detectable already by second generation interferometers |7j (see however ||). 

A generic prediction of unified theories (such as supergravity and superstring theories) 
is the existence of a gravitational multiplet which includes, beside the usual spin-2 graviton, 
scalar components. In a cosmological context, as well as the production of a primordial 
background of gravitational waves, it is therefore worth investigating the possible production 
of a relic background of scalar waves (for instance, dilatons |§). In particular, if the mass 
of the scalar field is small enough (m < 100 MeV), such a background would be present 
even today, and it could be accessible to direct experimental observation. It seems thus 
appropriate to estimate the sensitivity of gravitational antennas also to scalar waves. 

A relic background of scalar particles can interact with a gravitational antenna in two 
ways: either indirectly, through the geodesic coupling of the detector to the induced back- 
ground of scalar metric fluctuations, or even directly , through the effective scalar charge 
of the detector (such a direct coupling cannot be absorbed into the metric interactions, in an 
appropriate frame, if the scalar charge is non- universal). Up to now, the sensitivity to scalar 
waves has been mainly studied using the indirect coupling of the antennas to the scalar part 



of metric fluctuations, and considering massless scalar fields (see for instance fill]). Under 



those assumptions, the only difference with respect to the case of standard gravitational 



radiation is represented by the polarisation tensor of the scalar wave JT2[ . 

However, if the detector response tensor is characterised by a symmetric 3x3 trace- 
free tensor (such as the differential mode of an interferometer), the sensitivity to a scalar 
wave of momentum p and energy E generated by a massive scalar field is suppressed by a 
factor |T2j (p/E) 4 with respect to the sensitivity to an ordinary gravitational wave[|. This 
suppression is ineffective only for modes which are ultrarelativistic in the sensitivity band of 
the detector, m <^ E = (p 2 + m 2 Y^ 2 ~ p; this condition, however, occurs when the frequency 
corresponding to the mass of the scalar field is much smaller than the typical frequency of the 
detector, i.e. when m <C E ~ 10 — 10 3 Hz. The scalar field would be then associated to long- 
range interactions, and its coupling to the mass of the detector should be highly suppressed 
(with respect to the standard gravitational coupling), in order to agree with the existing 



tests of the equivalence principle and of macroscopic gravitational interactions [O. Hence, 



as far as interferometric antennas are concerned, the possible detection of a background of 



*This factor coincides with the suppression factor for the cross section of a resonant bar to longi- 



tudinal massive scalar modes [1C] 
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(massive or massless) scalar waves would seem to be strongly unfavoured with respect to 
the detection of a graviton background (in the literature, indeed, the possible detection of 
scalar waves is usually demanded to future resonant detectors of spherical type fT4||). 

Nevertheless, it is important to notice that, when the mass of the scalar field corresponds 
to a frequency which is in the detector sensitivity band, it is possible to obtain a resonant 
response also from the non-relativistic part of the scalar waves spectrum [Tj|. On the 
other hand, unlike a relativistic background of massless particles (like gravitons), the non- 
relativistic part of a relic background is not constrained by the nucleosynthesis bound, and 
could saturate the critical energy density required to close the Universe. As suggested in 
in ||15|| , if the non-relativistic branch of the spectrum is peaked at p ~ m, the polarisation 



suppression factor and the weakness of the scalar coupling could be compensated by the 
high relic density, and such a background of scalar waves could be a potential source for 
interferometric detectors. 

The aim of this paper is to discuss the possible detection of a relic stochastic background 
of massive scalar particles with the two LIGO interferometers (including the enhanced and 
advanced configurations). The paper is organised as follows. In Section |I| we recall the gen- 
eral expression for the optimised signal-to-noise ratio (SNR), obtained by cross-correlating 
two detectors, with respect to a stochastic background of massive scalar waves. In Section 
I fl] we apply this result to the case of the differential mode of an interferometer, taking into 
account both the geodesic coupling to the scalar part of the metric fluctuations, and the 
direct coupling to the scalar field. In Section [Iy| we discuss some examples, and we show in 
particular that a relic background of non-relativistic scalar particles, whose energy density 
provides a significant fraction of the critical energy density, can induce a non-negligible signal 
in the cross-correlation of the LIGO interferometers (in the enhanced and advanced config- 
urations), if the frequency corresponding to the mass of the scalar field is in the sensitivity 
band of the detectors. The main results of this paper are finally summarised in Section [V|. 



II. SIGNAL-TO-NOISE RATIO 

We will consider a stochastic background of massive scalar particles, described in terms 
of a scalar field <j)(x,t), and characterised by a dimensionless spectrum fl(p), 

tt{p) = — T — , Pc = — , (2.1) 

p c amp o7r 

where p is the momentum, p is the scalar field energy density, p c is the critical energy 
density, Hq the present value of the Hubble parameter, and Mp the Planck mass. We shall 
assume that the scalar field is coupled to the mass of the detector with gravitational strength 
(or weaker), and that the spectrum extends in momentum space from p to p\. The low 
frequency cut-off po may be zero (for growing spectra), or the present Hubble scale (for 
decreasing spectra), while p± is an high-frequency cut-off which depends on the details of 
the production mechanism. 

Our starting point is the expansion of the scalar field in the momentum space, 

(j)(t,x) = J™ dp J d 2 n{e 27rllE{p)t - ph ^<j){p,h) + h.c.} , (2.2) 
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where h is a unit vector specifying the propagation direction, p = hp is the momentum 
vector, and the energy E(p) of each mode is 



E(p) = ^p 2 + (m/2n) 2 (2.3) 

(we are adopting "unconventional" units h = 1, so that the proper frequency is simply / = 
E). The background of scalar waves is assumed to be isotropic, stationary and Gaussian || 
and satisfies the following stochastic conditions 

(<f>(p,h)) = o, 

{(pip, h)(f)*ip', h')) = 8{p - p')5 2 {h - ti)${p) , (2.4) 
where, using the explicit definition of 0(p), 

*M = i«L (2.5) 
87T d p E z [p) 

The scalar background induces on the output of the gravitational detector a strain h^t), 
proportional to the value of the scalar field <p(t,XD) at the detector position [f[2 



h^t) = J™ dp J dPh.F^h) {e 27Tl[E{p)t - ph -* D] <j){p,h) + h.c.} , (2.6) 

where is the position of the detector centre of mass, and F^ is the antenna pattern. In 
particular, 

F^h) =qe ab D ab (2.7) 

where e a b is the polarisation tensor of the scalar wave, D ab is the detector response tensor, 
and q is the effective coupling strength of the scalar field to the detector. The explicit form 
of F^ will be discussed in the next section. 

The optimal strategy to detect a stochastic background requires the cross-correlation 
between the output of (at least) two detectors 0, with uncorrected noises n»(i) , i = 1,2. 
Given the two outputs over a total observation time T, 

Si {t) = h^t) + m{t) , i = l,2, (2.8) 

one constructs a 'signal' S, 

rT/2 
-T/2 



r T/2 

S= / dtdt's 1 {t)s 2 (t')Q(t-t') , (2.9) 

J-T/2 



where Q(t — if) is a suitable filter function, usually chosen to optimise the signal-to- noise 
ratio: 

SNR=(S)/AS, (2.10) 

where AS 2 = (S 2 ) — (S) 2 is the variance of S. In our case, we can compute the mean value 
(S) by using the expansion ( |2.6| ) in momentum space, the statistical independence of the 



two noises (i.e. {ni(t)ri2(t')) = 0), and the fact that the noise and the strain are uncorrected 
(i.e. {ni{t)h l ^{t')) = 0). By assuming that the observation T is much larger than the typical 
intervals t — t' for which Q ^ 0, we obtain 

ff 2 ^ (too Q(E(p))Q(p) 7 (p) \ 
(5) = 5^ TRe U dp ^EHi) /' (2 - U) 



where 



and 



/oo 
dTe 27Ti{Eip)r Q(T) , (2.12) 



1{V) = 4^ / dVmph< " D1 '" D2) ^) F ^ ( 2 - 13 ) 

is the so-called overlap reduction function |3|], which depends on the relative orientation and 
location of the two detectors. In ( [2.131 ) the normalisation constant has been chosen so that - 
in the massless case - one obtains ^(p) = 1 for two coincident and coaligned interferometers. 

Switching to the frequency domain (E = f,dp/df = f/p), the mean value of S can be 
written as 

(s) A Re { r/(/--n)Q(/)»(v^h(v^W)| 

on 2 [Jo (j 2 — m 2 )j J 

where m = m/2ir, and 9 is the Heaviside step function. To compute the variance we will 
assume that for each detector the noise is much larger in magnitude than the strain induced 
by the scalar wave (i.e. rii(t) ^> h l ^{t)). One obtains J| 

AS 2 ^< S 2 >^ - jf dfP 1 (\f\)P 2 (\f\)\Q(f)\ 2 , (2.15) 
where Pi(\f\) is the one-sided noise power spectral density of the z-th detector, defined by 

MtfriW) = "? / dfe^-t )p i (|/|) . (2.16) 



£ J — oo 

Introducing the following positive semi-definite inner product in the frequency domain, 

(a,b) =Re{/ o °°a(/)&(/)Pi(/)P2(/)} , (2.17) 
the signal-to noise ratio can be written as 



(SNR) 2 = 2T 



H, 



2 



2 r 



(Q, A ) 



>V l(Q,Q)\ 

where 



(2.18) 



6(f - rh^i^P^W^^P^W) 

(P-m*)fP 1 (\f\)P 2 (\f\) ■ [ - J) 
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The above ratio is maximal if Q and A are proportional, i.e. Q = A A. With this optimal 
choice the SNR reads 

m 



SNR= f^-M V2TJ, (2.20) 



where 



T = dp ,, =4 - . (2.21) 

Jo Pi(y/p' 2 + m 2 ) P2^\Jv l + m 2 )(p 2 + m 2 ) 3 / 2 p 3 

For 7(p) = (df/dp)j(f) our result reduces to expression for the SNR already deduced 
in |U| (modulo a different normalisation of the overlap function). The scalar nature of the 
background is encoded into j(p), and will be discussed in the next Section. 



III. ANTENNA PATTERNS AND OVERLAP REDUCTION FUNCTION 

Given the spectrum, and the noise power spectral densities of the two detectors, the 
computation of the signal-to-noise ratio requires now the explicit expression of the overlap 
reduction function 7(p) for a pair of gravitational antennas. As already pointed out in 
the Introduction, we shall restrict our analysis to interferometric detectors. We will also 
consider the interaction of the scalar waves with the differential mode of the interferometer 



see e.g. |T2[), which is described by the following symmetric, trace- free tensor JT6 



Dab = ^{u a Ub- v a v b ), a, b = 1,2, 3, (3.1) 

where u a , v a are two unit vectors pointing in the directions of the arms of the interferometer. 
We will therefore neglect the interaction with the common mode, which is expected to be 
much more noisy. Note that in this paper we are not interested in distinguishing a scalar 
signal from a tensor one, but only in estimating the level of the signal eventually induced 
by a background of scalar waves. 

The computation of 7(p) requires knowledge of the antenna pattern ( |2.7| ), which de- 
scribes the induced strain and takes into account both the polarisation of the wave and the 
geometrical configuration of the detector (parametrised by D ab ). As far as the strain induced 



by a scalar wave is concerned, there are two possible contributions ||10| : one corresponds 
to the direct interaction of the detector with the scalar field, while the other is due to the 
indirect interaction with the scalar component of the metric fluctuations induced by the 
scalar field itself. 

Indeed, in a general scalar-tensor theory in which the matter fields are non-universally 
and non-minimally coupled to the scalar field (for instance, gravi-dilaton interactions in a 
superstring theory context jTTfl), the macroscopic bodies are characterised by a composition- 
dependent scalar charge (which cannot be eliminated by an appropriate choice of the frame, 
like the Jordan frame of conventional Brans-Dicke models), and their motion in a scalar- 



tensor background is in general non-geodesic [10|. Taking into account also the direct cou 



pling of a test mass to the gradients of the scalar background field, it follows that the standard 
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equation of geodesic deviation (which is the main equation for deriving the response of a 
gravitational detector) is generalised as follows fpLOTl : 

^ + R af3u *W£ a + ?£ a V a V0 = 0. (3.2) 

Here £ M is the spacelike vector connecting two nearby (non-geodesic) trajectories, and q is 
the scalar charge per unit of gravitational mass, representing the relative strength of scalar 
interaction with respect to ordinary tensor-type interactions |10| . 

The small, non-relativistic oscillations of a test mass, mechanically equivalent to the 
detector, are thus described by the equation: 

r = -e (R b00 a + qd b dy) . (3.3) 

The two contributions to the strain come from the gradients of the scalar field 0, and from 
the gradients of the scalar component if) of the metric fluctuations (induced by 0), covariantly 
represented by the Riemann tensor. The two fields <fi and ip are in principle different, but 
not independent, being related by a set of coupled differential equations (which are model- 
dependent). The antenna pattern F^(n), associated to the (indirect) Riemannian part of 
the strain, has been computed in \T2\ for both massless and massive scalar fields. We shall 
see that, for a traceless detector response tensor, the function F^{n) is also proportional to 
the antenna pattern F${n) associated to the direct, non-geodesic part of the strain. 

Introducing the transverse and longitudinal projectors with respect to the direction of 
propagation of the scalar wave, 



- ab 



(8 a b-h a h b ), L ab = n a n b , (3.4) 



the indirect, Riemannian contribution to the scalar pattern function becomes [12 

T ab + ^L ab J. (3.5) 

Using eq. (|3.3j ) and the mode expansion for the scalar field, the antenna pattern for the 
direct, non-geodesic coupling is: 

2 

F^n) = qD ab e ab ^) = qD ab ^ L ab . (3.6) 
Since T ab = 5 ab — L ab , and Tr D = 0, it follows that 

F^n) = -qF^n) = -q^F°(n), (3.7) 

where F® is the antenna pattern corresponding to a massless scalar wave The overlap 
reduction function of two interferometers, directly interacting through a charge g« with a 
scalar field, can thus be written as 

p \ 4 



lip) = Q1Q2 ( ) 7o(p) , (3-8) 
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where 70 (p) is the overlap reduction function for the geodesic interaction with a massless 
scalar field fll2|l . Using Eq. ( 2.21 ) and ( p.8|) , the signal-to-noise ratio is finally given by: 



SNR = qi q 2 

Note that this expression, with q^, — 1, is also valid to estimate the signal indirectly 
induced in the interferometers by the scalar metric fluctuations through their usual coupling 
to the Riemann tensor, provided f2(p) refers to the associated spectrum of scalar metric 
fluctuations. In the following we shall therefore use eq. ( |3.9| ) setting qi = q 2 = 1 if the 
dominant signal comes indirectly from the stochastic background of scalar metric fluctuations 
f^(p) induced by the scalar field, and setting instead g« < 1 (according to the experimental 
constraints) if the dominant signal is directly due to the stochastic background fi</,(p) of 
massive scalar waves. 



2T 



dp 



p 5 fi 2 (p)7o(p) 



Pi{Vp 2 + m ' 2 ) 



P 2 (Vp 2 + ™ 2 )(p 2 + m 2 ) n /2 



(3.9) 



IV. EXAMPLES 



We will now apply the results of the previous sections to estimate the signal induced by 
a stochastic background due to a massive scalar field on the two LIGO interferometers; we 
will consider the two detectors operating in the initial (I), enhanced (II) and advanced (III) 
configurations JT81. In particular, for the noise spectral density we will use the following 



analytical fits: 
• LIGO I @: 



P(f) = 



Pn = 10 



-16 




200Hz. 



(4.1) 



LIGO II [EDI: 




110Hz. 



(4.2) 



LIGO III |2 



P(f) 



5 



Jo, 



+ 2 + 2 



'£ 

Jo, 



P = 2.3 x lO^Hz" 1 , f s = 12Hz, f = 75Hz. 



(4.3) 
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In the above equations, f s is a seismic cut-off below which the noise spectral density is 
treated as infinity. We shall also assume that the non-relativistic modes in the spectrum 
are the dominant ones, and their energy density almost saturates the critical energy density, 
namely 

/ d\np n non - re \p) ~ 1. (4.4) 
Jo 

We shall analyse, in particular, two examples of spectra. 



A. Minimal dilaton background 

The first example is a stochastic background of massive relativistic dilatons, produced 
according to some models of early cosmological evolution based upon string theory || . To 
illustrate the difficulty in detecting such a background, we will consider here the dilaton 
spectrum obtained in the context of a "minimal" pre-big bang scenario and we shall assume 
that the dilaton mass is enough small so that the produced dilatons have not yet decayed 
into radiation with a rate T ~ m 3 /Mp. This implies T < Hq, i.e m 100 MeV. 

Even if the mass is negligible at the beginning of the radiation era, the proper momentum 
p = k/a(t) is red-shifted with respect to the rest mass because of the cosmological expansion, 
and all the modes in the dilaton spectrum tend to become non-relativistic. As a consequence, 
the present dilaton spectrum has in general at least three branches, corresponding to: i) 
relativistic modes, with p > m; ii) non-relativistic modes with p m < p < m, which became 
non-relativistic inside the horizon; in) non-relativistic modes with p < p m , which became 
non-relativistic outside the horizon. The separation scale p m corresponds to a mode that 
became non-relativistic, p ~ m, just at the time of horizon crossing, p ~ H, and is given 
by |^T| p m = pi^m/Hi) 1 / 2 , where Hi is the Hubble scale at the transition between the 



inflationary phase and the radiation phase, and pi ~ Hiai/a ~ (iJi/Mp) 1//2 10 11 Hz is the 
high frequency cut-off of the spectrum. 

The complete dilaton spectrum, at the present time t , can thus be parametrised as 
follows 



n(p) 



(it) 2 "?(*)(£)'> m<p<p u 

2 / 2 \ 1/2 / \ 5-1 



p m <p<m, (4.5) 



it, Ur ' p<p 



where fi 7 ~ 10 -4 is the present photon energy density (in critical units), if eq ~ 10 _55 Mp 
is the Hubble scale at the epoch of matter-radiation equilibrium, and 5 > is a slope 
parameter, depending on the kinematics of the pre-big-bang phase (5 = 3 for the low 
energy dilaton-driven phase, while 5 < 3 for modes crossing the horizon during the high 
curvature phase). 

Let us recall now that if the dilaton is strongly coupled to macroscopic matter with a 
composition-dependent charge q > 1, as suggested by loop corrections to the string effective 
action p3[, then the mass has to be large enough, m > 10~ 4 eV ~ 10 11 Hz, to be compatible 
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with known phenomenological bounds ||24j| . In that case, the mass is too far from the 
sensitivity range of a typical interferometric antenna, and no signal is detectable, since the 
spectral noises Pi are dominated by the mass, and tend to infinity. 

To discuss a possible detection, let us thus consider the alternative possibility of (almost) 



universal string loop corrections discussed in [Tj|, which is compatible with very light (but 



weakly coupled) dilatons, and let us assume the most favourable situation in which: i) 
the frequency corresponding to the dilaton mass lies within the sensitivity band of the 
interferometers; it) 5 < 1, and its value is enough small so that the spectrum is dominated 
by the non-relativistic branch; in) the spectrum saturates the critical density bound. In 
such a case, using the explicit form ( |4.5| ), it is easy to show that the main contribution to 
the integral ( [4.4|) comes from the low frequency band p < p m , and that the mass and the 
spectral index are related through the following equation: 

(H./Mp) 2 (m/f^) 1 / 2 (m/H^ 2 = 5. (4.6) 

With such a small mass, however, the dilatonic interaction is long-range and the direct 
coupling to matter is strongly suppressed (5 < 1 according to f24j]): the SNR turns out 
to be negligible, as shown by an explicit computation. Let's therefore make the further 
assumption that the dilaton background induces scalar metric fluctuations with a spectrum 
which has the same amplitude and the same frequency distribution as the dilaton spectrum 
( |4.5| ) (this is indeed what happens for relativistic dilatons, at the tree-level, thanks to the 
gravi-dilaton mixing of the string effective action ||). Such a background of scalar metric 
fluctuations would be gravitational coupled to the detectors with q — 1. 

We have computed the signal-to-noise ratio ( |3.9| ) obtained by cross-correlating the two 
LIGO detectors with q\ = q 2 = 1, using the background ( |4.5| ) with Hi fixed around the 
string scale Ms — H\ = 10 _1 Mp, and the slope S determined by the condition (^4.6|) . We 
have assumed T = 10 7 sec, and the overlap function 70 (p) has been computed using the 
coordinates of the two LIGO sites (see e.g. [p5] ). 



The results of this computation are summarised in Fig. 1, where the SNR is plotted as 
a function of the dilaton mass. The scalar signal would be detectable at a good confidence 
level by LIGO III, for a dilaton mass in the range 20 — 300 Hz. Unfortunately, however, 
such values of m compatible with detection would imply for the metric fluctuation spectrum 
a maximum intensity Q ~ 1 around a peak scale p m ~ 10 -10 — 10~ 9 Hz. It follows from eq. 
( [4.5D that Q > 10 -2 at 10 -8 Hz, in sharp contradiction with the bound on metric fluctuations 
obtained from pulsar-timing data p6| (which requires Q < 10~ 8 at 10~ 8 Hz). Imposing the 
pulsar-timing constraint on the induced spectrum of scalar metric fluctuations it turns out 
that the signal is not detectable. 



B. Non-relativistic scalar dark matter 

The difficulties in detecting the minimal dilaton background discussed in the previous 
section are mainly due to the fact that the spectrum (f4.5|) is characterised by a peak located 
V = Pmi much lower than the mass scale p = m. However, if the peak of the 
spectrum is not too far from the mass scale, the signal could be detected even through 
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the direct interaction of the scalar background with the interferometers, with an effective 
coupling compatible with the existing phenomenological bounds. 

To illustrate this possibility we shall therefore consider a stochastic background of scalar 
massive particles whose energy density represents a significant fraction of the critical energy 
density. We will assume that the spectrum is dominated by non-relativistic modes, with a 
peak for p ~ m. For the sake of simplicity we use the following single-power law for the 
spectrum: 

n(p) = 5 tt s (-Y , p<m, 5>0, (4.7) 

where the parameter Q s is the total energy density (in units of the critical energy density) of 
the spectrum, integrated over all modes. Since we are interested in a mass for the scalar field 
corresponding to a frequency in the range 10 — 10 3 Hz, we shall assume that the scalar charge 
q is consequently suppressed: taking into account all existing phenomenological bounds [fJI| . 
it turns out that in this frequency range the allowed coupling can be parametrised as 

logg 2 < - 7 + log(m/10 Hz), 10 Hz < m < 1 kHz (4.8) 

for universal scalar interactions, and 

logg 2 < - 8 + log(m/10 Hz), 10 Hz < m < 1 kHz (4.9) 

for composition-dependent interactions. 

We have computed the signal-to-noise ratio, obtained by cross-correlating the two LIGO 
detectors, for the spectrum ( |4.7|) assuming a nearly critical energy density (Q s = 1) and 
the maximal value of q allowed by the bound ( (4.8[ ). Fig. 2 shows the contour plots of the 
SNR for T = 10 7 function of log 10 (m) and of the spectral index S, for the initial (I), 

enhanced (II) and advanced (III) LIGO configurations. We have considered the case 8 < 3, 
typical of string cosmology models || . Already with the enhanced LIGO configuration, the 
signal produced by this class of backgrounds may have a significant signal-to-noise ratio 
in the relevant mass range. Similar results can be obtained with a composition dependent 
scalar charge, provided the values of SNR are rescaled by a factor 10 _1 , according to eq. 

Fig 3 shows the values of q 2 Q s corresponding to SNR=1 for the enhanced and advanced 
LIGO configurations, as a function of the scalar mass, using the spectrum ( |4.7[ ) with 5 = 1 
(which is a typical value for the slope of the tensor perturbations spectrum predicted in 
quintessential inflation models p7l). The interesting result is that, in the absence of any 
signal, an upper limit on q 2 Q s would imply an indirect bound on the scalar charge stronger 
than the direct bounds existing at present f24|, provided the expected energy density of the 
scalar background corresponds to a non-negligible fraction of critical energy density. 

V. CONCLUSIONS 

In this paper we have discussed the possible detection of a relic background of mas- 
sive, non-relativistic scalar particles, through the cross-correlation of two interferometric 
gravitational antennas. 
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Using previous results [|T(],|T2|,|T^] , we have computed the signal-to-noise ratio for a stochas- 
tic background of massive scalar waves, taking into account both the scalar and gravitational 
charge of the detector, i.e. including in the differential mode of the interferometer both the 
direct coupling to the scalar field spectrum and the coupling to the induced spectrum of 
scalar metric fluctuations. 

We have computed the signal-to-noise ratio obtained by cross-correlating the two LIGO 
interferometers in the initial, enhanced and advanced configurations, for two particular ex- 
amples: relic dilatons produced in the context of string cosmology models, and a generic 
scalar component of dark matter. 

If the total energy density of the scalar background is dominated by non-relativistic 
modes and amounts to a significant fraction of the critical energy density, we have found 
that the induced signal is certainly non-negligible for the LIGO enhanced and advanced 
configurations, provided both the mass and the peak of the spectrum lie within the sensitivity 
band of the detectors, i.e. from 10 Hz to 1 k Hz. These results hold even if the direct coupling 
to the detector is strongly suppressed, according the existing bounds on long range scalar 
interactions. 

It is worth pointing out that the differential mode of an interferometer is particularly 
unfavoured for the detection of non-relativistic scalar waves: since the detector response 
tensor is traceless, the corresponding pattern function, for non-relativistic modes, is sup- 
pressed by a factor (p/E) 2 ~ (p/m) 2 <C 1. This is the reason why, to obtain a detectable 
signal, we have restricted our discussion to spectra peaked around p = m. However, for 
different detectors (for instance, resonant spheres) the detector response tensor with respect 
a scalar wave is not traceless, and the suppression factor is absent. It would be interesting 
to investigate the detectability of a more general class of stochastic backgrounds of massive 
scalar particles (i.e. where the peak and the mass are not in the same frequency range) 
through resonant spherical detectors. 

The main conclusion of our work, confirming the results of previous studies, is that a 
stochastic background of light but non-relativistic scalar particles, which provides a sig- 
nificant fraction of critical energy density, could induce a non-negligible signal in future 
enhanced interferometric detectors (in competition with the signal of a possible stochastic 
background of gravitational waves), provided the frequency corresponding to the mass of the 
scalar field lies within the detector sensitivity band. We believe that this possibility should 
be taken into account in the data analysis development of second-generation interferometric 
detectors. 
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FIG. 1. The signal-to-noise ratio ( ]3.£| j for the cross-correlation of the two LIGO interferometers, 
coupled with maximal strength (qi = q2 = 1) to the scalar spectrum (4.t), which is normalised 
so as to saturate the critical density at p = p m . The integration time is T = 10 7 sec. The signal 
is possibly non-negligible only if the mass is in the sensitivity band of the gravitational detectors. 
It must be noted, however, that the extrapolation of the spectrum to much lower scales is not 
compatible with the bounds obtained from pulsar timing data. 
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FIG. 2. Contour-plots in the plane {m, 6} for the SNR of eq. 



, computed for the scalar 
spectrum (4. / ), and relative to the cross-correlation of a pair of LIGO interferometers in the initial 
(I), enhanced (II) and advanced (III) conhgurations. We have used T = 10 sec, Q s = 1, and the 
maximal allowed value of q according to eq. ( [4,^ j. The region inside the contour corresponds to a 
value of SNR bigger than the value on the border. 
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FIG. 3. Allowed values of q 2 Q s , as a function of the mass of the scalar held, and corresponding 
to SNR = 1, computed for the scalar spectrum ( |4, /[ ) and for the two LIGO detectors in the 
enhanced and advanced conhgurations. For comparison, we have also plotted the maximal value 
of q 2 according to eq. ( |4.gp (solid line). We have used a total observation time T = W 7 sec, and a 
spectral index 6 = 1. 
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